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ABSTRACT 

Aims. A near-infrared study of the main ionizing star of the ultracompact HII region G23.96+0. 15 (IRAS 18317-0757) is presented, 
along with a re-evaluation of the distance to this source, and a re-assessment of H- and K-band classification diagnostics for O dwarfs 
Methods. We have obtained near-IR VUT/ISAAC imaging and spectroscopy of G23.96+0. 15, plus archival imaging from 
UKIRT/UFTI. A spectroscopic analysis was carried out using a non-UTE model atmosphere code. 

Results. A quantitative H- and K-band classification scheme for O dwarfs is provided, from which we establish an 07.5V spectral 
subtype for the central star of G23.96+0. 15. We estimate an effective temperature of Teff ~38 kK from a spectral analysis. 
Conclusions. A spectroscopic distance of 2.5 kpc is obtained for G23.96+0. 15, substantially lower than the kinematic distance of 
4.7 kpc, in common with recent studies of other Milky way H n regions. Such discrepancies would be alleviated if sources are unre- 
solved binaries or clusters. 
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>■ 1. Introduction 
O 

P\| The formation of high mass stars remains an unresolved prob- 

lO • lem in contemporary astrophysics (Zinnecker & Yorke 2007). 

CO ■ Despite considerable theoretical progress in recent years, this 
topic remains observationally challenging. O stars form deeply 

^-H embedded within their natal cocoons, presumably within inter- 

OO mediate to high mass star clusters, only to visually emerge after 
~0.5 Myr (Prescott et al. 2007), a a substantial fraction of their 
2-10 Myr main-sequence lifetimes. 

• I— I Radio continuum surveys of ultracompact H n (UCHII) re- 
gions betray the presence of O stars through the effect of 
Jh their Lyman continuum photons to free-free (thermal) emission. 

. . ■ Typically, these regions can not be observed at wavelengths 
shorter than the mid-infrared, such that the ionizing stars can 
only be studied indirectly through the circumstellar dust and 
gas (e.g. Peeters et al. 2002; Martm-Hernandez et al. 2002). 
However, in a few instances our sight line to the central source 
is sufficiently clear of circumstellar dust that the O star respon- 
sible for the H ii region can be directly observed at near-infrared 
wavelengths (Watson & Hanson 1997; Hanson et al. 2002; Bik 
et al. 2005). 

Over the past decade, the advent of efficient near-IR spectro- 
graphs at large ground-based telescopes has permitted medium 
resolution (R ~ 5000) spectroscopy of template Milky Way O 
stars in the H- and K-bands (Hanson et al. 2005b), which has 
been extended to the 04-5 star responsible for the G29.96- 
0.02 UCHII region (IRAS 8434-0242, Hanson et al. 2005a). 
Spectroscopic analysis of near-IR hydrogen and helium lines 
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of O stars agrees closely with optical diagnostics in most cases 
(Repolust et al. 2005). 

Thus far, G29.96-0.02 represents the sole example of an 
UCHII region whose ionizing star has been analysed based upon 
medium spectral resolution. The significance of such objects 
is that, in principal, they may serve as calibrators for the so- 
called 'inverse problem', in which the ionizing stars of em- 
bedded UCHII regions may be obtained from analysis of mid- 
infrared fine-structure nebular lines (e.g. Martm-Hernandez et 
al. 2002; Morisset et al. 2004). 

Fortunately, other UCHII regions are also accessible to near- 
IR spectroscopy. One such region, G23.96-h0.15 (IRAS 18317- 
0757) is the focus of the present study, for which Hanson et al. 
(2002) estimated a spectral type of 07-8 from low resolution 
K-band spectroscopy and Kim & Koo (2001) have inferred a 
mid-O spectral type from radio continuum observations for an 
adopted distance of 6 kpc. Hunter et al. (2004) have identified a 
cluster of embedded massive stars from millimetre observations 
ofG23.96+0.15. 

In the present study, we present new near-IR imaging and 
spectroscopy of G23. 96-1-0. 15 permitting the subtype of the ion- 
izing star to be refined from H-band and K-band hydrogen and 
helium line ratios, plus a spectroscopic distance for comparison 
with kinematic results. 



2. Observations 

H-band and K-band spectroscopy of G23. 96-1-0. 15 was observed 
with the ISAAC instrument (Moorwood et al. 1998) mounted at 
the Very Large Telescope between 22 ApriI-2 May 2006. High 
spatial resolution imaging was drawn from acquisition ISAAC 
datasets, which were supplemented by archival UKIRT Fast 
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Fig. 1. (left) 2x2 arcmin colour composite JHK image cen- 
tered upon G23.96+0.15 from 2MASS; (right) 10x10 arcsec 
~2.18;um acquisition image of G23.96+0.15 from ISAAC, in 
which the brightest sources are marked with the ISAAC slit over- 
laid. Source #1 has position a = 18:34:25.25, S = -07:54:45.5 
(J2000), with the locations of the 6 cm peak from Wood & 
Churchwell (1989, circle) and 21cm peak from Kim & Koo 
(2001, square) also indicated. North is up and east is to the left. 



Trace Imager (UFTI, Roche et al. 2002) K-band imaging, plus 
H-band imaging from Hanson et al. (2002) and 2MASS datasets. 

2.1. Near-IR imaging 

Acquisition ISAAC images of G23.96h-0.15 were obtained using 
the 1024x1024 Hawaii Rockwell array (0.148 arcsec/pix), with a 
combination of both the 2. 17jum and 2. 19jum narrow-band filters 
during excellent seeing conditions of 0.3-0.5 arcsec in April- 
May 2006. 

In Fig.[T|we present a 2x2 arcmin JHK composite image cen- 
tered upon G23.96-h0.15 from 2MASS plus the central 10x10 
arcsec of the ~2.18//m ISAAC dataset. The ISAAC image re- 
veals a bright source (#1), located at a = 18:34:25.25, 5 = - 
07:54:45.5 (J2000.0). This source Ues -2.7 arcsec west of the 
6 cm peak reported by Wood & Churchwell (1989), though 3.5 
acsec east of the 21 cm peak of Kim & Koo (2001) for this irreg- 
ular UCHII region. Fainter near-IR sources lie 1.5 arcsec to the 
north of #1 and 4.2 arcsec to the south, which we shall refer to as 
sources #2 and 3, respectively. Source #3 is itself resolved into 
NW and SE components, separated by 0.7 arcsec. Photometry of 
these sources was obtained from archival UKIRT UFTI images 
from 19 Jul 2000 obtained with the 1024x1024 Hawaii aiTay 
(0.09 arcsec/pix) and K-band (K98) filter during seeing condi- 
tions of ~0.6 arcsec, using a zero-point obtained from 10 nearby 
2MASS field stars, providing an accuracy of +0.04 mag. The J-, 
H- and K- band magnitudes for the combined sources #1 and #2 
are 13.81, 11.46 and 9.94 mag, respectively (2MASS 18342523- 
0754455). 

We have also inspected H-band images from the Steward 
2.3m IRCam guider images of Hanson et al. (2002) using a 
128x128 Hawaii array (0.5 arcsec/pix), obtained in seeing con- 
ditions of ~1.5 arcsec, in which sources #1 and #2 are again 
blended, from which H(#3) - H(#l+#2) ~ 1.5 mag. Independent 
estimates of H(#3NW)-H(#2) ~ 0.3 mag and H(#3NW)-H(#1) 
~ 2.3 mag have been obtained from our high spatial resolu- 
tion ISAAC 1.71/zm long-slit spectroscopy (see next section). 
These sources were well aligned, so slit losses should minimal. 
Photometric properties are presented in Table [1] taking into ac- 
count uncertainties resulting from the ISAAC and IRCam obser- 
vations. 



Table 1. Near-IR properties of G23. 96-1-0. 15 sources from 
UKIRT/UFTI (K-band) and Steward 2.3m/IRCam imaging (H- 
band) plus VLT/ISAAC acquisition images and spectroscopy, 
from which equivalenth widths (Wa in A) were also measured 
for#l. 



Source K H-K W^CHeii) IV^(Hei) lV,i(Hi) lV,i(Heii) Spect. 

mag mag 1.692yum 1.700yum 2.165yum 2.189yum Type 

1043 L2 020 04 235 094 07.5V 

±0.04 ±01 ±0.08 ±011 ±017 ±012 

#2 12.00 1.7 

±0.04 ±0.2 

#3NW 12.10 1.9 

±0.04 ±03 

#3SE 13.7 -02 

±0.2 ±03 



2.2. Near-IR spectroscopy 

Long-slit spectroscopy of G23.96+0.15 was obtained with 
ISAAC in April-May 2006 at a position angle of 7.7 degrees 
west of north, in order to include sources #1, #2 and #3NW at 
three medium resolution, 0.775A/pix) grating settings centered 
at 1.71, 2.09 and 2.20//m (recall Fig.[TJ. These observations, ob- 
tained using a 0.6 arcsec slit width, were taken during excellent 
seeing conditions (0.3-0.5 arcsec) at low airmass (1.02-1.08) to- 
gether with solar-type telluric analogues. 

Six individual exposures, comprising three AB pairs, were 
obtained for each grating setting, with wavelength solutions 
achieved from comparison XeAr arc datasets. From these, 
the observations covered 1.671-1.751jum, 2.029-2. 155jum and 
2.140-2.265jum at spectral resolutions of 3.8A, 6. OA and 6. OA 
respectively, as measured from arc lines. 

Telluric correction was achieved from a single AB pair of 
spectroscopic datasets of early-G dwarfs observed at a similar 
airmass to G23.96+0.15, corrected for their spectral features us- 
ing high resolution observations of the Sun which were adjusted 
to both the radial velocity and spectral resolution of the template 
stars. 

The 2.09jum setup suffered from low-level variable struc- 
ture which was accentuated upon flat-fielding, so that only the 
two other settings were flat-fielded. Consequently, the contin- 
uum S/N achieved was ~100 for the 1.71/vm setting and 150 for 
the 2.20;um setting for #1, but only 60 for the 2.09;um setting. 
The continuum S/N for sources #2 and #3NW was, at best, no 
greater than 40-50. 

3. The ionizing star of G23.96+0.15 

3.1. Near-IR spectral classification 

Hanson et al. (1996) developed a classification scheme for O 
stars from low resolution K-band spectroscopy. Classification 
of early to mid-O stars was achieved from the presence of C iv 
2.08//m and Nin 2.11yum emission lines, although it proved to 
be problematic to distinguish between late O and early B stars. 
This technique was extended to the H-band by Hanson et al. 
(1998) who highlighted the diagnostic role of the hydrogen Br 
11 1.681yum, He II 1.692//m and He: 1.700/im lines. 

In Fig. |2] we present H-band and K-band spectroscopy of 
G23.96+0.15 #1 together with MK classification O dwarfs from 
the high quality, medium resolution atlas of Hanson et al. 
(2005b). From visual inspection, a classification of 07-8 is es- 
timated, in agreement with the low resolution K-band spectrum 
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Fig. 2. H-band (left) and K-band (right) VLT/ISAAC specti-oscopy of G23.96+0.15 #1 together with MK template 06-9.5V stars 
from Hanson et al. (2005b), together with a spectral fit to #1 (dotted line, see text) 



of #1 from Hanson et al. (2002). However, we have also inves- 
tigated the possibility of using the Hen 1.692yum/Hei l.lQQfim 
and He ji 2.189;um/Br'x ratios for quantitative classification, in- 
dependent of the metal lines in the K-band (see also Lenorzer et 
al. 2004; Repolust et al. 2005). 

We have measured the equivalent widths of these lines for 
MK classification O stars from Hanson et al. (2005b), whose 
absorption line ratios are presented in Figl3] (excluding Bry 
emission line supergiants). From this it is clear that the Hen 
1.692/im/Hei 1.700/im ratio provides an excellent classifica- 
tion diagnostic for O stars, especially dwarfs, while the Hen 
2. 1 89//m/Bry ratio provides an additional constraint for late-type 
O dwarfs, providing the spectral resolution is sufficient to ex- 
clude any nebular contamination. Polynomial fits to equivalent 
width ratios of dwarf stars from Hanson et al. (2005b) are pre- 
sented as a guide, excluding HD 37468 from the fit for which 
negligible He n absorption is observed. The H-band ratio inher- 
ently represents a superior diagnostic, although the K-band ratio 
has two advantages for late-type O dwarfs, namely a lower inter- 
stellar extinction plus 2.189/im is an intrinsically stronger Hen 
line than 1.692yum. 

Unfortunately, the moderate S/N of sources #2 and #3 
prevented identification of spectral feaures beyond Bry nebu- 
lar emission. We present H and K-band photometry of these 
sources in Table [1] together with absorption line measurements 
of G23.96+0. 15 #1 . From the lower S/N H-band ratio we infer a 
dwarf spectral type of 08^,' j while the K-band diagnostic indi- 
cates O7.5+0.5. Therefore, we propose an 07.5V classification 
for source #1. 



3.2. Stellar temperature of G23.96+0. 15 #1 

The current spectroscopic calibration of O stars would result in 
an effective temperature of reff=36kK for an 07.5 star (Martins 
et al. 2005). Instead, we have obtained a direct estimate of the 
stellar temperature from a comparison between the near-IR hy- 



drogen and helium lines and synthetic spectra obtained with the 
CMFGEN model atmosphere code (Hillier & Miller 1998). 

CMFGEN solves the radiative transfer equation in the co- 
moving frame, under the additional constraint of statistical equi- 
librium. The temperature structure is determined by radiative 
equilibrium. Since CMFGEN does not solve the momentum 
equation, a density or velocity structure is required. For the su- 
personic part, the velocity is parameterized with a classical /3- 
type law, with an exponent of /3-1 adopted. This is connected 
to a hydrostatic density structure at depth, such that the velocity 
and velocity gradient match at the interface. The subsonic veloc- 
ity structure is set by a corresponding fully line-blanketed plane- 
parallel TLUSTY model (v.200, see Lanz & Hubeny 2003). The 
atomic model is similar to that adopted in Hillier et al. (2003), 
including ions from H, He, C, N, O, Ne, Si, S, Ar, Ca and Fe. 

We have assumed a depth-independent Doppler profile for 
all lines when solving for the atmospheric structure in the co- 
moving frame, while in the final calculation of the emergent 
spectrum in the observer's frame, we have adopted a uniform 
turbulence of 50 km s ' . Incoherent electron scattering and Stark 
broadening for hydrogen and helium lines are adopted. Finally, 
we convolve our synthetic spectrum with a rotational broadening 
profile. 

For an adopted terminal wind velocity of 2000 kms"', sur- 
face gravity of log g - 4 and abundance ratio of He/H=0. 1 by 
number, we varied the stellar radius and (non-clumped) mass- 
loss rate until a reasonable match to the Hen 1.692yum, Hei 
1.700;um, Bry and He n 2.189;um was achieved. An acceptable fit 
was achieved for Teff = 38 +lkK, logMAMoyr"') = -6.3+0.2 
and vsin/ ~ 220 + 30 kms"', as shown in Fig.|2] Nebular con- 
tamination is significant for He i 1 .700/im and high members of 
the hydrogen Brackett series. For an adopted mass of 30 Mq, a 
radius of 9.2 Rq results, which is representative of normal Milky 
Way mid-O dwarfs (e.g. Repolust et al. 2005). We defer any fur- 
ther discussion of the physical properties of G23. 96-1-0. 15 #1 un- 
til the next section where its distance is considered. 
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Fig. 3. (Top) Equivalent width ratio of Hen 1 .692//m/He i 
1.700yum versus subtype for O dwarfs (circles), giants (squares) 
and supergiants (triangles) from Hanson et al. (2005a) and 
G23.96+0.15 #1 (filled circle) together with a polynomial fit to 
dwarfs as a guide (solid line); (bottom) as above for ratio of He ii 
2.189/im/Bry 

Table 2. Physical properties of G23.96+0.15 #1 from our spec- 
tral analysis, based upon either the spectroscopic (top) or kine- 
matic (bottom) distances. Mass estimates are approximate since 
we adopt logg=4.0. The predicted Lyman continuum ionizing 
flux (Qo) is shown together with lower limits obtained from the 
21 cm radio observations of Kim & Koo (2001, ™') 
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4. Distance and extinction to G23.96+0.15 

Distances to Galactic Hii regions are typically obtained from 
kinematic methods, using an adopted rotation curve (Brand & 
Blitz 1993) and Solar galactocentric distance (8 kpc, Reid 1993). 
Following this approach the near (far) kinematic distance to 
G23.96-h0.15 is 4.7 kpc (lO.Okpc) based upon the LSR velocity 
of 79.3 kms"' from Wink et al. (1983), suggesting a distance of 



4.2 kpc from the Galactic Centre for G23.96-h0.15. LSR veloci- 
ties from NH3 (Churchwell et al. 1990), CS (Plume et al. 1992) 
and H76q' recombination lines observations (Kim & Koo 2001) 
agree with Wink et al. (1983) to within +1 kms"'. 

One would expect Mk = -3.9 mag for G23.96-h0.15 #1 
(07.5V) from the Conti et al. (2008) absolute visual magnitude- 
spectral type calibration, together with the intrinsic colours of 
Martins & Plez (2006). We may correct its observed K-band 
magnitude for extinction using A K - 1.82^q 23£'h-k (Indebetouw 
et al. 2005) and (H - K)o = -0.10, from which Ak = 2.33 + 0.3 
mag, suggesting a distance modulus of 12.0^^^, i.e. a distance of 
2.5+0.4kpc. For a representative scatter of +0.5 mag in the ab- 
solute magnitude calibration of O stars, a spectroscopic distance 
in the range 2-3.2 kpc would be implied, placing G23.96-h0.15 
at a distance of 5.8+0.5 kpc from the Galactic Centre. 

Martm-Hernandez et al. (2002) have previously estimated 
Ak - 2.0 mag for G23.96+0.15 based upon HI recombina- 
tion lines observed in Infrared Space Observatory (ISO) spec- 
troscopy. We independently confirm this result from the ob- 
served ISO Brff/Br/J ratio, although the ISO pointing was off- 
set by >10 arcsec from the source peak, and so measured fluxes 
need to be treated with caution (Peeters et al. 2002). 

In common with other studies of H 11 regions, we find a spec- 
troscopic distance to G23.96+0.15 that is substantially lower 
than kinematic distances (e.g. Blum et al. 2001; Figueredo et al. 
2008). The two distances could be reconciled if #1 were substan- 
tially more luminous than typical mid-O stars, with Mk = -5.5 
mag, although this would require the star to be extremely mas- 
sive if it is single. This is highlighted in Table |2l where physical 
properties of G23.96+0.15 #1 for each of the alternate distances 
are presented. Of course, #1 could be an unresolved equal mass 
binary (or compact star cluster). Such a source would appear 
~0.75 (0.3) mag brighter in the K-band than a single 07.5 star, 
owing to the contribution of its companion(sfl A binary/cluster 
scenario would help to alleviate the spectroscopic and kinematic 
distance discrepancy, although in the cluster case, stellar absorp- 
tion lines from the O star would be diluted by the continuum 
light from lower mass cluster members, which does not appear to 
be the case for G23.96+0.15 #1 (recall Fig.|2]i. Similar remarks 
apply for G29.96-0.02, for which a spectroscopic distance of 
~3 kpc would be expected for a single 04 dwarf using near-IR 
photometry from Pratap et al. (1999), compared with a kinematic 
distance of 7.4kpc (Fish et al. 2003). 

Kim & Koo (2001) estimated an ionizing output of 10"*** ** 
ph s ' for G23.96+0.15A from their 21 cm radio continuum ob- 
servations, which is equivalent to an 07 dwarf (Conti et al. 
2008), albeit based upon an assumed distance of 6 kpc. If one 
was to adjust the distance to 4.7 (2.5) kpc, a reduced Lyman con- 
tinuum output would be obtained, representative of an 08 (09.5) 
star (Table |2]|. The measured 21 cm flux, ~0.75 dex lower than 
our indirect estimate, is merely a lower limit since it may be par- 
tially optically thick to bremstrahlung radiation and a fraction 
of the ionizing photons are likely to be absorbed by dust (either 
within the H 11 region or circumstellar cocoon). 

The near-IR colours of sources #2 and #3NW are consistent 
with being physically associated with G23.96+0.15 #1 (#3SE 
is likely a foreground source). For the spectroscopic distance 
of 2.5 kpc, they are probable early B dwarfs, at projected dis- 
tances of 0.02 pc (#2) and 0.05 pc (#3NW) from#l. In common 
with indirect dust and gas diagnostics for other UCHII regions 
(Okamoto et al. 2003), we find evidence that G23.96+0. 15 likely 



' On the basis of a Salpeter-like Initial Mass Function for high mass 
stars 
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hosts multiple massive stars (see also Hunter et al. 2004). For 
the kinematic (near) distance of 4.7 kpc, sources #2 and #3NW 
would possess absolute magnitudes typical of late O stars, and 
so would also most likely be unresolved sub-clusters within the 
UCHII region, at projected distances of 0.035-0.1 pc from #1. 

5. Conclusions 

We show that G23.96+0.15 provides a second UCHII region 
whose ionizing star (#1) is accessible to high S/N, medium res- 
olution near-IR spectroscopy, complementing the higher ion- 
ization G29.96-0.02 source. We provide a quantitative near- 
infrared classification scheme for O stars, which has potential 
application for other visibly obscured O dwarfs in the Milky 
Way, from which a near-IR 07.5V classification is obtained for 
G23.96-)-0.15 #1. We obtain a stellar temperature of reff=38 kK 
for G23.96+0. 15 #1 and infer a spectroscopic distance of 2.5 kpc 
for a single star origin. This is substantially smaller than the 
(near) kinematic distance of 4.7 kpc, in common with other ob- 
scured H II regions (e.g. Blum et al. 2001 ; Figueredo et al. 2008), 
although an unresolved binary or cluster nature would help to al- 
leviate this discrepancy. Two fainter sources are also likely early- 
type members of the UCHII region. The availability of spec- 
tral properties for G23. 96-1-0. 15 should help empirically address 
the 'inverse problem' of obtaining the ionizing O stars for other 
UCHII regions based solely upon mid-IR nebular characteristics. 
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